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I. INTRODUCTION 

The theoretical analysis of pinning phenomena of elas- 
tic objects in random potentials is an important physical 
problem with a great impact on many fields of research 
[ll, [1] . Elastic objects in disordered media subject to a 
constant driving force are meanwhile well understood [sl 
at zero temperature and the influence of finite tempera- 
ture has been studied as well [4]. 

In recent years, also the problem of ac-driven elastic 
interfaces in disordered systems has gained experimen- 
tal interest 0-01 • experiments, considerable attention 
is devoted to the behaviour of the ac-susceptibility of 
ferroelectric thin films or ultrathin ferromagnetic multi- 
layers, which is believed to be related to the motion of 
domain walls in these systems. So, an interest to under- 
stand ac-driven domain wall motion in disordered sys- 
tems emerged. First results towards an explanation of 
the hysteretic behaviour of the magnetisation [1, |^ have 
been achieved and the observed frequency regimes for 
the response of ferroic domain wall motion are under- 
stood phenomenologically jioj . An important contribu- 
tion comes from the velocity hysteresis of moving domain 
walls fllj , but the characterising parameters have not yet 
been worked out. Generally, with regard to an analytic 
description of the domain wall motion characteristics, re- 
sults are scarce. Yet, mainly numerical results have been 
obtained to describe qualitative features of such systems, 
like the hysteresis or double- hysteresis of the velocity [iTI - 

In many cases, for a first quantitative analysis of com- 
plicated non-linear problems one uses a perturbation ex- 
pansion. However, perturbative approaches are some- 
times hampered by mathematical subtleties, like non- 
analyticities or singular perturbation theory (cf. [T^). 
or by physical obstacles such as non-perturbative excita- 
tions or strong coupling. The difficulty with perturbation 



theory as a tool for the analysis of pinned elastic objects 
has its own interesting history. 

Until the beginning of the eighties, the lower critical 
dimension di of the random field Ising model has been 
the subject of a long-lasting debate. Dimensional reduc- 
tion predicted that the lower critical dimension equals 
d/ = 3, whereas domain- wall arguments lead to the 
conclusion, that = 2. Eventually, in 1984 a final 
decision could be made and dimensional reduction was 



proven to fail [1 
been found later 
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The reason for the failure has 
1^ to be connected with the ex- 
istence of many metastable states for the domain walls 
separating different regimes in a multidomain configura- 
tion. The plethora of metastable states arises from the 
dominance of the disorder over the domain wall elastic- 
ity on length scales above the so-called Larkin length Lp 
[20| for sample dimensions d < 5. A perturbative iter- 
ation to find the energy minimum will not necessarily 
yield the correct extremal state [111, [l^]. Put in more 
mathematical terms, the formal perturbative treatment 
of the domain walls assumes an analytic disorder correla- 
tor. However, a functional renormalisation group (FRG) 
treatment shows [l8l |. that any initially analytic disorder 
correlator develops a cusp-singularity at a finite length 
scale, which is the Larkin length Lp. 

This insight has important consequences for the prob- 
lem of an interface in a disordered environment exposed 
to a constant driving force h. At zero temperature and 
for small external force ft., the interface adjusts its con- 
figuration to balance the driving and the disorder, but 
remains pinned and does not move on large time scales. 
If h is tuned to exceed a critical threshold /ip, then after 
transience has relaxed, the interface slides with a mean 
velocity v that behaves as w ~ {h — hpY for h \ /ip. The 
system undergoes a non-equilibrium phase transition, the 
so-called depinning transition, for which v plays the role 
of an order parameter. The critical behaviour close to the 
depinning transition has been investigated in a number of 
works [23-_30j. Furthermore, the influence of finite tem- 
perature on this transition has been considered in several 
articles [1, |3l|, [s^l ■ Even the scaling behaviour for small 
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frequencies lo of an ac-drive on approaching the critical 
point {h,uj) = (/ip,0) has been worked out 

The depinning transition can, however, not be ac- 
counted for by perturbatively expanding the disorder. 
Therefore, perturbation theory is not applicable for small 
constant driving forces. 

To achieve deeper quantitative insight into the ac- 
dynamics of interfaces, perturbation theory seems to be 
the only feasible method. Clearly, in the vicinity of the 
critical point — (/ip,0) perturbation theory is ex- 

pected to yield erroneous results, if not properly com- 
bined with an FRG treatment. Yet, FRG equations for 
the similar problem of a constant driving force, have been 
obtained by the construction of a perturbative series and 
a subsequent e-expansion. So, understanding the pertur- 
bation theory is the first necessary starting point. As has 
been said, pure perturbative calculations are restricted 
to parameter regimes far away from the neighbourhood 
of the depinning transition. Thus, for sufficiently large 
frequencies and driving field amplitudes, from the phys- 
ical point of view there appears to be no contraindica- 
tion against a perturbative procedure. However, as will 
become clear in this article, for systems subject to pe- 
riodic driving forces perturbation theory in the disorder 
strength gives non-regular contributions to the velocity 
corrections for an internal interface dimension D < A. 
With the attribute non-regular we refer to expressions 
that grow unboundedly in time. Such unbounded con- 
tributions certainly do not reflect the true physical be- 
haviour, but their origin deserves careful investigation. 
Far away from the critical point, this behaviour of per- 
turbation theory cannot be related to the assumption of 
an analytic disorder correlator. Quite the contrary, work- 
ing with a cusped disorder correlator brings additional 
difficulties due to the delta functions in its derivatives. 

In section |TT1 we introduce the equation of motion for 
ac-driven elastic manifolds in disordered media, and in- 
vestigate its perturbative expansion. Although our model 
is taylored to describe elastic manifolds, like interfaces 
between two immiscible fluids or domain walls in fer- 
roic systems, we believe that our analysis also covers a 
wide range of models for other interesting problems, e.g. 
charge density waves [ssj or flux lines in type-II supercon- 
ductors jSJ-jSa]- After analysing the flrst non- vanishing 
order, we derive the diagrammatic expansion to account 
for higher orders. This can be used to argue, that per- 
turbation theory works for D > A. The failure of pertur- 
bation theory for D < 4 is then analysed and explained. 
The well-known suitability of perturbation theory for es- 
timates of the velocity of domain walls driven by a con- 
stant force, far in the sliding regime, does not contradict 
our statements for the ac-driving. We are going to take 
a look at this as well. The technically involved and more 
mathematical treatments are taken out of the main text 
and given in the appendices [B] and [Cj 

The failure of perturbation theory in all physically in- 
teresting cases D < 4 underlines the importance of the 
mean-field approach, which is the second central subject 



of this article. In section Hill we explore the correspond- 
ing mean-field equation of motion and its perturbative 
expansion. After some illustration of the qualitative be- 
haviour of the full solution, we prove, that the pertur- 
bative corrections are regular, i.e. they remain bounded 
in all orders. The bulk part of this inductive proof, the 
induction step, is outsourced to appendix [El Further, we 
show that for large enough driving field amplitudes, suf- 
ficiently strong elastic coupling and high frequencies, the 
perturbative results agree very well with the numerics for 
the full mean field equation of motion. We conclude our 
considerations of the mean-field problem by a numerical 
analysis of the lowest non-vanishing perturbative order, 
which yields the decay law of the Fourier coefficients with 
the driving field strength. 

A table of frequently appearing symbols is provided in 
appendix [F] 



II. THE PERTURBATIVE TREATMENT OF 
INTERFACES IN DISORDERED MEDIA 

A. The model 

Our analysis models interfaces and domain walls that 
are thin such that they can be described by elastic D- 
dimensional manifolds, embedded in a I? -|- 1 dimen- 
sional space. The manifold itself is parametrised by a D- 
dimensional set x of coordinates and its position in space 
is given by z{x,t). We confine ourselves to the study of 
the zero-temperature case, i.e. we do not take thermal 
noise into account. Moreover, our model assumes small 
gradients and does not allow for overhangs. We expose 
the interface to a periodic driving force 



h{t) = h ■ cos ujt. 



(1) 



Then, the overdamped dynamics of elastic interfaces can 
be described by the equation of motion 



j-^dtz{x, t) = rVlz + h{t) + u ■ g{x, z), 



(2) 



which has already been introduced in earlier works (stI - 
[39| . In eq. ((21), F and 7 are the stiffness and the in- 
verse mobility of the domain wall. For simplicity, we set 
7 = 1 in the following. The function g{x^ z) describes the 
quenched disorder, taken to be Gaufiian with the corre- 
lators given by 



{9{x,z)) 
{9{x,z)g{x',z')) 





S'^ix 



c')A{z-z'), 



(3) 
(4) 



where (...) denotes the average over disorder. The disor- 
der correlator in z-direction is taken symmetric around 
and decays exponentially on a length scale £. Further, we 
demand A(0) = 1, as the strength of the disorder shall 
be measured by u. To be definite, we choose 



A(z — z') — exp 



(5) 
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in case we need a precise formula. This choice corre- 
sponds to the case of an elastic manifold exposed to ran- 
dom field disorder [l^. Throughout the whole paper we 
assume weak disorder. This means, that pinning forces 
are weak and the interface is pinned at the fluctuations 
of the impurity concentration, and not at single pinning 
centres. A more precise definition can be found e.g. in 
(40| . For weak disorder, the random forces have to ac- 
cumulate to overcome the elasticity. On small length 
scales, elastic forces dominate and the interface is essen- 
tially fiat. By comparing the elastic and the disorder 
term in ([2]) one can estimate the length scale Lp, called 
the Larkin length, at which the two competing effects are 
of the same order. The result is 



re 



(6) 



Thus, the elastic term dominates on all length scales for 
£) > 4. 

Finally, we specify the initial configuration for the 
equation of motion ([2]) to be a fiat wall z{x, t = 0) = 0. 



B. Perturbation theory - first order 

The equation of motion ([2|) can only be solved via an 
expansion in the disorder strength u. We are going to de- 
rive the perturbation expansion directly for the equation 
of motion, since this appears simpler. There is, however, 
another approach via functional integrals which came in 
useful for the functional renormalisation group calcula- 
tions in the case of a constant driving force. The reader 
who is more familiar with this technique may find a brief 
treatment in appendix]^ which reveals the connection to 
our methodology. 

The expansion is naturally performed around the solu- 
tion for the problem without disorder, i.e. where u — 0. 
In this case, we have a flat wall following the driv- 
ing field: Z{t) — {h/uj)smujt. Thus, we decompose 
z{x,t) into the disorder-free solution and a correction, 
i.e. z{x,t) = Z(t) + (^{x,t). The equation of motion for 
the disorder correction C{x,t) is easily derived from eq. 



m 



{dt-TVl)ax,t)=u-g{x,Z + C). 



(7) 



The Green function for the differential operator on the 
left hand side is the well-known heat kernel 



idt-TVl)Gix,t) 
G{x,t) 



S^{x)6{t) 

m 



d^fc 



Akx-rk^t 



The fc-integral has to be cut off at some scale A, corre- 
sponding to the inverse smallest length scale in the sys- 
tem. To set up the perturbation series, we expand the 
correction in the disorder strength 



and the disorder force around the non-disordered solution 
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{x,Z + 0^T.^29(.^,Z)^ 



k=0 



(10) 



Thus, we obtain an equation for the first order correction: 

oo 

Ci(a;,t)= J <^"x' j dt'G{x~x',t^t')g{x',Z{t')). (11) 



Obviously, the disorder average vanishes. The disorder 
average for the second order contribution is given by 



(C2) (t) ^J<^t^J di2 A'[Z(ii) - Z(i2)]x 


-rfc2(ti-t2) 



d^yfc 



(12) 



The second order correction to the velocity follows 
straightforwardly 



{V2)it)^ / dt' A'[Z{t)-Z{t')] 



(27r)^ 



(13) 



To get a first impression on how this expression behaves 
for large t, we split off the Fourier-O-mode: 



A'[Z{t) - Z{t')] = Fo{Lot)/£ + p{t, t') 



(14) 



n=0 



F^iiot) = ^ ( _ ) . sin(2n + l)u:t. (15) 



Here, p{t,t') is a well-behaved oscillation around in t' . 
Since p{t,t') yields a bounded contribution to {V2) (t), 
we consider only Fo{t) to find out, how {V2) increases 
asymptotically in time. The Fourier coefficients Kn can 
be determined analytically. They are diminished when 
their argument increases or approaches zero and they re- 
main bounded. Integration over t' yields 



Foiiot) Sd 



dk 



(27r)^ J r 




,D-3 



1 - e 



-k^rt 



t— Fojiot) 



where 



1-e-P 



(16) 
(17) 

(18) 



ax,t) 



(j)^ For X 00 the integral ajj converges for D < 2 and 
diverges logarithmically for D — 2. Thus, the asymptotic 
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FIG. 1: Plot of the first non-vanishing perturbative correction 
{«2) {t) to the disorder average {v{x, t)) for different interface 
dimensions. For the plot we used h — 1 and the units are 
chosen such that uj = i = 1. 



behaviour of the first perturbative correction in time is 
given by 




(19) 



where coit) is some bounded function. Fig. [T] shows 
the plots of {V2) {t) for D = 1, 2, 3. Obviously, there is a 
problem of the perturbation expansion in low dimensions. 
In the following, we are going to see that perturbation 
theory does not work even for D < 4. 



C. Higher order graphical expansion 

Higher orders of the perturbation expansion are best 
expressed diagramatically. To deduce the diagrammatic 
rules, all one has to do is plugging ([9]) into (fTO|) . rear- 
ranging the sum in powers of u and inserting this into 
([7]). The diagrammatic rules that emerge are fairly sim- 
ple: to express the perturbative correction of order n, we 
draw all rooted trees with n vertices and add a stem. Up 
to the fourth order, this tree graph expansion is given by 



-• • 



< 



-• • • 



(20) 
(21) 

(22) 



-• • • • 



< 




(23) 



Every vertex represents a disorder insertion 
f;(™)(Z(t))/m!, where m counts the number of out- 
going branches (away from the root). Every line 
corresponds to an integral operator, the kernel being the 
propagator G{x,t). To get the graphical expansion for 
the velocities, just remove the first line. The disorder 
average can be carried out using Wick's theorem, since 
our disorder is assumed to be Gaufiian. An immediate 
consequence is, that after averaging over disorder, only 
graphs with an even number of vertices survive. In the 
following, we are going to consider the perturbation 
expansion for the disorder averaged velocity v{t). The 
pairing for the disorder average shall be denoted by a 
dashed line. An example graph from the 4*''-order is 

•■ (24) 



Using this graphical expansion, in appendix [B] we take 
a look at the general behaviour of the diagrammatic con- 
tributions to all orders for t ^ 00 and argue that all 
graphs remain bounded in case D > A. 



D. The failure of perturbation theory for D < 4 

For D < 4: the perturbation expansion is not as well- 
behaved as for D > 4. In this section, we show this 
and give an explanation why perturbation theory is ill- 
behaved in low dimensions. 

Though the disorder averaged graphical structure can 
become complicated, one especially simple graph has the 
same structure in all orders: the one for which all ver- 
tices are connected directly to the root. In the equations 
(|20ll23p . we have drawn those graphs at the very first 
place. Let us call them bushes. The general bush graph 
contribution to the velocity correction of order 2p thus 
corresponds to the following diagram 



2p) 




(25) 



where the dotted line is a placeholder for other vertices 
that we have not drawn. Up to combinatorical factors, 
the general (disorder averaged) bush i?2p that occurs in 
the 2p-th order perturbative correction to the velocity v, 
reads 



{B2p)it) (X 



dtidt2 A[Zi - Z2] 







t 

X J dt' A^^P-^'^[Z - Z'' 


0-1 



(27r)^ 



-rfe^(2t-ti-t2) 



^-rfe^(t-t') 



(26) 
(27) 



5 



To work out the asymptotic envelope of {B2p) (t) for large 
t we only need to take a look at the expression Ti . Since 
A is an entirely positive function, and the fc-integral is 
certainly non-negative, we can replace A by its maximum 
A(0) to get an upper bound, and if we replace A by 
the minimal value that is taken A[2h/{Lu£)], we get a 
lower bound for Ti. In both cases, A[Z{ti) — Z{t2)] is 
replaced by a constant, so that the integration over ti and 
t2 can be done easily. Consequently, up to a constant, 
the asymptotic envelope of Ti is given by 



£2 



Ti{t) 



u 



d^fc 



1 



-D 



(28) 



Here, the constants t and d are time scales, given by 





4-D 





u 







Ll/r 



(29) 
(30) 



and the function A£){x) has been introduced for nota- 
tional convenience 



Ad{x) = 



{2n) 



dpp 



D-5 



1-e-P 



(31) 



The function Ajj (x) is increasing but bounded for D < 4 
and grows logarithmically as a; ^ oo for D — A. Thus Ti 
grows monotonically in t for any D < A. 

Actually, our statement about the asymptotics in eq. 
([28|) is more robust than our simple argument suggests, 
and in fact does not rely on the positivity of A. A more 
detailed calculation, presented in appendix [Cl shows that 



u 



Ti = 



t 

T 

u^A^ 



KD 



2p 



Pd 



uA 2 

t 



(32) 



and 



y2^2(p-l) 



uA 2 








u^A" 


PD 





ID 



-, Lot 



(33) 



where all of the functions kd, ku, Pq, fjj and pu are 
bounded in t for D < 4. The first term of Ti shows, that 
for D < A the whole expression grows for t — > oo without 
any bound. This is, because by definition (cf. equation 



(jClOl) ). K_D a Ad, where Ad is given by ((3T|) . The other 
tree graphs that appear in the graphical representation of 
{v2p) and that we have not analysed here, exhibit similar 
behaviour. Cancellations among diagrams do not occur. 
To have a little more evidence, that D — A really enters 
as an upper critical dimension, we have perturbatively 
investigated the interface's width in appendix|D] In sum- 
mary, perturbation theory is ill behaved for D < A. The 
reason for this shall be discussed in the following. 

The bushes {B2p) that we have considered so far are 
part of an expansion of the disorder averaged velocity 



v(t) — h ■ cosujt + Vdis{t) 



(34) 



in the disorder strength u. The dimensionless ratios, in 
which u occurs in that expansion are 



and 



uA 2 



(35) 



Since in the stationary state the interface is expected to 
follow the driving with frequency w, its Fourier represen- 
tation has to take on the form 



cosnujt+ (36) 




s'mnujt. (37) 



Because of the ratio t/r, an expansion in the disorder 

u brings powers of in every order, since r depends 
on u (cf. equation ([291) '). The remaining question is the 
meaning of r. Since r appears as a time scale for the time 
dependence of the Fourier coefficients, which physically 
should approach a constant value in the stationary state 
t — ?► oo, the most natural interpretation is the transience. 
Keeping in mind, that we start with a flat wall, we have 
to expect several kinds of transience effects. As we have 
seen, there are only two time scales in question: r and t?. 
As can be concluded from their definitions (cf. equations 
and (PU)) ). they obey the same structure: an intrin- 
sic length scale to the power 2 divided by T. The time 
scale T involves the Larkin length Lp (cf. equation ©), 
which measures the competition between disorder and 
elasticity: the interface is flat on length scales L < Lp. 
This indicates that, up to some dimensionless prefactor, 
T describes the time during which correlated interface 
segments of extension Lp adopt to their local disorder 
environment, i.e. the roughening time of the interface. 
For D > A, the interface is flat on all length scales, thus 
there is neither roughening nor a Larkin length, hence t 
is meaningless and cannot occur. This agrees with our 
observation. For D > A, there is no disorder-dependent 
time scale any more, which could bring powers of time in 
the perturbative corrections, therefore they remain finite 
a,s t oo. However, also for D > 4 the disorder leads to 
a typical deviation of every point of the interface from the 
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mean position. The built-up of this typical deviation to- 
wards its steady-state value happens on the time-scale i), 
in agreement with the observation for mean-field theory 
(cf. section UlI Cp . Thus, both time scales can naturally 
be interpreted as the life times of two different transience 
effects. 



E. Interfaces subject to a constant driving force 



In the introduction, we already pointed out that per- 
turbation theory in connection with interfaces driven by a 
constant force cannot properly account for the existence 
of the depinning transition and therefore gives misleading 
results for D < 4. But far above the depinning thresh- 
old, i.e. for h ^ hp, the interface slides and its velocity 
can be estimated perturbatively. The dynamical corre- 
lation length Ly — Ti/v [l^l is then small compared to 
the Larkin length Lp and thus working with an analytic 
disorder correlator and expanding the disorder in its mo- 
ments works. Of course, also for constant driving forces 
one has to start with a definite initial condition, which is 
usually a flat wall. Thus, there will be transience effects 
for dc-driven interfaces as well [4l], \^ . In this section we 
take a short glance to understand the difference between 
ac and dc-driving. Our special interest is devoted to the 
time scales that determine the duration of transience ef- 
fects. 

The equation of motion for the elastic interface expe- 
riencing a constant driving force 



dtz{x, t) — rV^z + /i + u • g{x, z), 



(38) 



has the disorder-free solution {u = 0) Z{t) = ht. The 
perturbation expansion is essentially the same as in sec- 
tion lll B1 just the non-disordcd solution around which we 
expand is different. 

Actually, there is a problem with the decomposition 
z — Z+( here, since the sliding velocity v is different from 
h, hence C {v — h)t is not a small quantity (compared to 
£) for large t and the Taylor expansion (ITOl) of the disorder 
is questionable. Since here we shall not be interested in 
large times t > £/{h — v) but only want to determine 
the time scale of the transience (occuring at small t ^ 
i/{h — v)), this problem is safely ignored. 

The first non-vanishing correction to the velocity is 



found to be (cf. eq. ([T3)) ) 



(^^2) (t)= / dt' A'[Z{t)-Z{t')] 



-rfc^(t-t') 



(39) 



1 r d^fc 



h J (27r)^ 



_j.2pj_,x2t2 fTik^ ht 

e — (f 



2h 



(40) 



Sp 








(p i 


htV 


\2ht' 





(41) 



where we have introduced the function 



ip{a, 6) = 1 - \/^a • e"' • [erf(a + b) - erf(a)] (42) 

for convenience. The time-scales on which transience ef- 
fects disappear are obviously given by 



£ A~^ 
Tdc — — = — and d — — 

h T r 



^ (43) 



and are manifestly disorder-independent. Lh = Tt/h de- 
notes the correlation length of an interface moving with a 
velocity h. Of course, asymptotically the interface moves 
with a velocity v < h, but at the very beginning, when 
we start off with a flat wall, its velocity is indeed given 
by h (cf. dSHl)). 

It is not surprising, that the time scale of the initial 
roughening is different for dc and ac driving. Either prob- 
lem involves completely different physical processes to be 
responsible for transience. In the case of an ac-driving, 
the system undergoes a process of adaption of its configu- 
ration to the local disorder, such that a stable stationary 
oscillation is possible, and during which higher Fourier 
modes build up. For dc-driving, the system starts to 
move with a velocity of h, which then rapidly decreases, 
and roughens since segments of the interface are pinned 
and remain at rest until they are pulled forward by the 
neighbouring segments through the elastic coupling. The 
time for this process mainly depends on the velocity of 
the interface, not on the strength of the disorder. 



III. MEAN FIELD THEORY 

To extend the study of ac-driven interfaces beyond nu- 
merics, perturbation theory seems unavoidable. Since, 
as we have seen in the preceeding section, perturbation 
theory only works for D > 4, we turn here to an inter- 
esting limiting case, formally corresponding to D = cxi: 
the mean-field equation of motion. 

Mean field calculations have been performed before for 
the problem of interfaces and charge density waves sub- 
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jec t to a constant driving in a number of articles [2J, |43 
|4y|. Special emphasis has been put on the depinning 
transition and its critical properties. 

The perturbation expansion for dc-driven interfaces 
has been investigated by Koplik and Levine [s^, who 
also emphasised on the mean-field problem. They al- 
ready mentioned the same problematic graphs in their 
expansion that we will encounter below, but they did not 
provide a proof for the fact, that the unbounded terms 
cancel to all orders, independent of the driving. 



A. The mean-field equation of motion 

The mean field equation corresponding to our original 
equation of motion ([2|) is obtained via the replacement 
of the elastic term by a uniform long-range coupling (cf. 
e.g 



4J]). To do this, we have to formulate the model 1^ 
on a lattice in a;-direction, i.e. the coordinates that pa- 
rameterise the interface itself are discretised. The lattice 
Laplacian reads 



^ z{xi + acd) + zjxi - aed) - 2z(xi) 

D N 



XI 51 -^^3^ bi^Jd) ~z{Xi)\, 



d=l3d = l 



(44) 



(45) 



where a denotes the lattice constant. To get the mean 
field theory, Jij has to be replaced by a uniform coupling 
but such that the sum over all couplings "Y^- Jij remains 
the same. Hence, we choose 



J, 



MF 



1 



(46) 



Now, the disorder has to be discretised as well, which is 
achieved if we replace the delta function in the correlator 
dl by S'^^x^ -Xj) S^ja-^/"^ (cf. (H). The result- 
ing equation of motion should be independent of x, just 
the lattice constant a and the dimension enter because 
the disorder scales with a factor a^^/^. Finally, for the 
mean-field equation of motion, we obtain 



dtz 



[{z)~z] + h{t)+r^-g{z), 



(47) 



where c — T /a? and rj — ujaP^'^ . The disorder remains 
Gaufiian with 



(.9(^)) = 
{g{z)g{z')) = A{z^z'). 



(48) 
(49) 



The function A(z — z') is as before, so we shall choose 
again ([5]) whenever we need an explicit expression for 
calculations. 

The physical picture of the mean field equation of mo- 
tion is a system of distinct particles, moving in certain 




FIG. 2: Numerical solution of (|47|) for different driving field 
strengths and c = 1.0, rj = 2.5. For the simulation, t and z 
are measured in units such that u = £ = 1. 



realisations of the disorder. All of them are harmoni- 
cally coupled to their common mean, i.e. the elastic cou- 
pling between neighbouring wall segments TV^z is now 
replaced by a uniform coupling c • [(z) — z] to the disor- 
der averaged position (z), which in turn is determined 
self-consistently by the single realisations. 

Apart from the correlation length £ of the disorder, 
there is another important length scale in the system. In 
the absence of any driving force (i.e. ft. = 0), we can easily 
determine the mean deviation of the coordinate z of a 
special realisation from the disorder averaged position 
(z). For ft = we expect z = 0, at least in the steady 
state and ((Tf|) straightforwardly leads to 



(((^ 



(50) 



So, rj/c measures the order of the average distance from 
the common mean. 

In what follows, the disorder averaged velocity v = (z) 
will be denoted by the symbol v. 



B. Qualitative behaviour and numerical results 

To get an idea about how the system, corresponding 
to the equation of motion with an ac-driving (cf. (HTl) ) 



h{t) = ft • cos utt 



(51) 



behaves, we implemented a numerical approach. The 
disorder is modelled by concatenated straight lines, the 
values of the junction points are chosen randomly from 
a bounded interval. The correlator has been checked to 
be perfectly in agreement with ([5]). 

Before discussing the numerical trajectories, we note a 
first property of the equation of motion (j47p . It contains 
a symmetry of the (disorder averaged) system, namely 
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FIG. 3: Numerical solution of equation (|47p for h = 6.0, 
c = 1.0 and -q — 2.5 for different frequencies, t and z being 
measured in units such that ujq = £ = 1- 

that all disorder averaged quantities are invariant under 
the transformation h ^ —h and z — > — z, which implies 
V — > —v. We have hereby fixed the initial condition to 
be z(0) = for all realisations. If one chooses another 
initial condition, its sign has to be inverted as well, of 
course. In the steady state, i.e. for t (as is the 

time-scale on which transience effects are diminished, see 
below), the trajectory must therefore obey the symmetry 
h — >■ —h, V — >■ —V. This symmetry is obviously reflected 
in the numerical solutions (see fig. [5]). 

An interesting consequence of this symmetry is, that 
the even Fourier coefficients of the solution v{t) (which is 
periodic with period ^-k jio) vanish. Once the steady state 
is reached, the symmetry requires v(i) — —v(t + tt/uj). 
For the even Fourier modes this means 

2-71 

C2N = J dt v{t)e'^^'^' 



= Jdt ?;(t)e'2^"* + J Atv{t+ ^) e'2^"* = 0. 



(52) 

The typical picture of a w-/i-plot is that of a single hys- 
teresis for h <^ T] and a double hysteresis ioi h ^ rj. In 
an intermediate range, we find a single hysteresis with 
a cusped endpoint. The qualitative shape of the solu- 
tion trajectories, examples of which are shown in fig. [21 
agrees with numerical results [ill, [13 , that have been ob- 
tained as solutions for © in the case of finite interfaces 
with periodic boundary conditions. Moreover, as the fre- 
quency is sent to zero w — > 0, the hysteretic trajectory 
approaches the depinning curve for an adiabatic change 
of the driving field. This is shown in fig. [31 

In the following, we want to give a qualitative discus- 
sion of the hystereses in the case of small elasticity c, 
more precisely our discussion assumes r]/c^ £. 

Weak fields h <^r]: For each disorder realisation, the 




FIG. 4: Numerical solution of equation (|47|) for different elas- 
tic constants and h = 40.0, 77 = 10.0. The units of t and 2 are 
chosen such that to = £ = 1. 

configuration is allowed to deviate from the mean by an 
order oirj/c (cf. (l50l) ). which by our assumption is large 
compared to the disorder correlation length I. So, except 
for rare events, in the case of weak driving fields, the typ- 
ical system in a certain disorder configuration remains in 
a potential well of the disorder. Starting at some large 
time to (m the steady state) for which /i(io) — 0, we ex- 
pect a certain realisation to be located close to a zero 
point zq of a falling edge of the disorder force 3(20) — 0, 
since this corresponds to a stable configuration. As the 
field grows, the system starts to move in the direction of 
growing z, where the disorder force competes with the 
driving. Because in the vicinity of the potential mini- 
mum, the disorder force g{z) behaves approximately lin- 
ear in z, the acceleration is approximately zero and the 
velocity almost constant. This changes when the driv- 
ing is about to reaching its maximum. The slower the 
growth of the driving, the smaller the velocity. At the 
maximum, the velocity equals zero, as the driving and 
the restitutional disorder force compensate. For decreas- 
ing h{t), the restitution force wins and pushes the system 
back in the direction of the potential minimum. Hence, 
the velocity v turns negative short-time after the field 
has reached its maximum and is still positive. Once the 
stable position zq is reached again, the same starts in the 
negative direction. 

Certainly, the restitutional disorder force need not con- 
tinuously grow with z, but may exhibit bumps or similar 
noisy structure, but those details average out when tak- 
ing the mean over all disorder configurations. 

Strong fields h ^ rj: In the case of strong driving 
amplitudes, we encounter the situation of a double hys- 
teresis. Again, starting at a time t^ with /i(to) = for 
to 3> c~^, we may assume the system to be located close 
to a zero g{zQ) — oi a. falling edge of the disorder force 
field. As h{t) grows, we first have the same situation 
as in the case of weak driving: The disorder acts resti- 
tutionally and thus keeps the velocity small and leads 
to a small slope dv/dh. Once the field is of the order 
rj, the system is no longer locked into a potential well. 



9 



but a cross-over to sliding behaviour sets in. On fur- 
ther increasing h, the system finaUy arrives at a slope 
s = dv/dh, which depends mainly on i] and h. The 
larger h, the closer s approaches the value of 1. After the 
field reaches its maximum, the velocity decreases with 
the field, the slope being smaller than s, if h is not too 
large. This is explained as follows. The mean deviation 
^((z) — z)^) that we have estimated to be of the order 
r/'^/c^ for h ~ (cf. eq. ((5(I)) ). under driving also depends 
on the difference between v{t) and h{t) (as can be seen 
from (|47p ). At the beginning of the crossover to sliding, 
when h ^ T], this difference is large and therefore also the 
mean deviation is large (especially for small c). When h 
increases, the particles in disorder realisations that are 
behind the mean are strongly accelerated. On the other 
hand, those that are ahead of the mean position are not 
so much decelerated, because for each disorder realisa- 
tion, being ahead of or behind the mean position rapidly 
changes and a motion backwards (against the driving) is 
suppressed. Thus the reduction of the mean deviation 
towards some asymptotic value gives a stronger slope for 
rising fields during the crossover. When the field h{t) 
reduces, the mean deviation shrinks, h(t) approaching v, 
and V drops less rapidly than it has risen because still 
the mean deviation prefers a diminution of the difference 
v{t) — h[t). This approximately linear reduction of the 
velocity remains, until the field is weaker than the typi- 
cal disorder force, when the system is again trapped in a 
potential well. Since on rising edges of the disorder force, 
driving and disorder point in the same direction, the sys- 
tem will rarely sit there (it moves away very fast). The 
velocity becomes negative before /i = 0, since the system 
slides down the falling edge {dg/dz < 0) of the disorder 
force. At ft, = everything starts again in the negative 
direction. An example for fairly large field amplitudes is 
shown in fig. 21 

So far, our discussion has emphasised on small c. The 
effect of larger c is to couple the configuration z{t) of 
every realisation strongly to the mean {z{t)). This wipes 
out the effect of disorder. Thus for larger c the double 
hysteresis winds around a straight line, connecting the 
extremal velocities. This can be seen in fig. |4l 



C. Mean- field perturbation tlieory 

1. The diagrammatic expansion 

Since the mean-field equation of motion (1471) cannot be 
solved exactly, we attempt an expansion in the disorder 
strength rj. Therefore, as before, we decompose z = Z + 
C, where Z{t) = {h/uj) smujt is the solution of the non- 



disordered problem (ry = 0) around which we expand, 
and 

oo oo 

C-ECfc'?'. (C)-E(0fc'7'- (53) 

fc=l k=l 

is the perturbative correction. Still, we have the equa- 
tions for Cfc depending on (C)fc, which is also unknown. 
This eventually leads us to a set of two coupled equations 

{dt+c)C = c{0+Ti-g{Z + (54) 
dt{C.)=V-{9{Z + 0), (55) 

that we can solve iteratively for every order of the per- 
turbation series, if we expand 



(56) 



Tl = 



If one is interested to keep small orders, this expansion 
of the disorder can only work if C <C because i is the 
typical scale on which g{z) changes. We will come back 
to that point later in section ITlI C 3| when discussing the 
validity of perturbation theory. For the moment, we just 
do it. 

The propagator corresponding to the left hand side of 
Eq. dMl) reads 



G{t) = Q{t) ■ e" 



(57) 



Using this propagator, we can formally write down the 
solution and express it order by order in a power series 
in rj. Up to the second order, the solutions are 



(C)i(0 = o, 

Ci(i)= / dhe-<'-'^^g{Z{h)), 



(58) 
(59) 

(C)2(t)= / dh I dt2e-<'^-'^^A'[Z{h)-Zit2)], (60) 





t 

C2{t) = J dtie-^(*-*^) [c (C)2 ih) + g'izih)) ■ Ci(ti)] . 



(61) 

Since we assume Gaufiian disorder, the disorder averaged 
corrections {C)n vanish for odd n. We use a diagrammatic 
representation to depict the nested perturbation expan- 
sion. For the interesting quantities (C)^, the first two 
non- vanishing orders are given by: 
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(C)4= 3 




(62) 



•> • • 



r 



The diagrammatic rules are fairly similar to those in 
section III CI we draw all rooted trees with k vertices, 
and add a stem. Each vertex corresponds to a factor 
(7(™)(Z(i))/m!, where m counts the number of outgoing 
lines (away from the root). The line between two vertices 
represents a propagator G{t). Then Wick's theorem is 
applied to carry out the disorder average. Each two ver- 
tices, that are grouped together for the average, will be 
connected by a dashed line. Finally, there is one new fea- 
ture, that we did not come along in section III CI Every 
straight line which, upon removing it, makes the whole 
graph falling apart into two subgraphs, has to be replaced 
by a curly line. A curly line symbolises the propagator 
of ([55|) . which is just a Heaviside function &{t). Those 
graphs that contain an internal curly line are exactly the 
one-particle reducible (IPR) diagrams. 



2. Regularity of the perturbative series 

The perturbation expansion leaves some questions, 
that have to be addressed. It is not immediately obvious. 



J 



that taking the disorder average of (I5T|) gives the result in 
((60)) . i.e. (C)2 {t) = {C2{t))- However, a short calculation, 
using integration by parts reveals this relation to hold. 

Another, much deeper problem is related to the dia- 
grams involving a curly line in their interior. Due to the 
curly line, they grow linearly in time. Already in sec- 
tion III CI and appendix [B] we have mentioned that there 
are trees that contain lines the assigned momentum of 
which equals 0. Here, for the mean-field problem these 
lines are found as the troublesome curly lines: they con- 
nect a subtree with internal GauiBian pairing. Koplik 
and Levine [39| explicitly checked for a time independent 
driving h{t) = ft, up to sixth order, that the problematic 
terms of the IPR diagrams mutually cancel. We give a 
very general version of this proof, that holds for any time- 
dependence of the driving field h(t) and covers all pertur- 
bative orders. To illustrate, how this works, we present 
the calculation for the fourth order here. The somewhat 
technical induction step, which extends our argument to 
all orders is given in appendix [E) For simplicity, we work 
with the velocity diagrams, that are obtained by just re- 
moving the curly line from the root. 



t t-2 







dti e-^(*"*i)A"[Z(t) - Z(ti)] / dh / dhe~'^'^~'''^^'[Z{t2)-Z{h 



dti e 





tl t2 



-c(t-tO(_A"[Z(t)-Z(ti)]) j dt2 j dt, e-^(*^-*^)A'[Z(t2) - ^(is)] 





t tl 



S^jdhj dt2 e-^(*-*^)A"[Z(t) - Z{t2)] j e-=(*^-*^'A'[Z(ti) - Z^)] 



I 



(63) 

(64) 
(65) 
(66) 



The modification of the second diagram to express it as the sum of the first and 5* is merely integration by parts 
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for the integral over ti. The term S now corresponds 
to the sum of the two diagrams. It is easy to see, that 
S remains bounded for large times. Every time integral 
carries an exponential damping term. Basically, we have 
thereby established, that the perturbation series exists 
and is well-behaved in the sense, that there are no terms 
that lead to an overall unbounded growth in time. 



3. Validity of perturbation theory 

Still, the question is open, whether one may assume ( 
to be small compared to £. This was a requirement for 
the Taylor expansion ([56l) to be valid. If c is large, any 
particle moving in a particular realisation of a disorder 
potential is strongly bound to the disorder averaged po- 
sition. This prevents it from exploring the own disorder 
environment and thus large c effectively scale down 77. 
All realisations stay close to the disorder averaged po- 
sition, the mean deviation being approximated by 77/c. 
A problem now occurs, if the disorder averaged position 
deviates strongly from the 77 = solution. For ft. 3> 77 this 
can only happen during those periods, where h(t) takes 
on small values. The time that has to elapse until every 
system has adopted to its own disorder realisation, and 
hence the time until the system can be pinned, is 
(see below). For perturbation theory to work, this time 
must be large compared to the length of the period dur- 
ing which h < r], which we roughly estimate as ri/[ioh). 
This gives us a second condition for the applicability of 
perturbation theory: h/rj 3> c/uj. 

In summary, the conditions for perturbation theory to 
hold are the following. The driving force amplitude h 
has to be large compared to 7], h/rj ^ max{c/w, 1} to 
make the series expansion work and to guarantee that 
the disorder averaged solution stays close to the 77 = 
trajectory (around which we expand). Moreover, c must 
be large (c 3> "n/ 1) to ensure proximity of each realisation 
to the disorder average. 

The diagrammatic prescription yields up to terms of 
order 0{Tf) 



){t) cos ujt + if j di'e-'=(*-*')A'[Z(i)-Z(<')]. (67) 




A direct comparison of the numerical solution of (|67p and 
the full equation of motion ()47|) , shown in fig. [SJ confirms 
an excellent agreement. 



D. Perturbative harmonic expansion 

For an ac driving force, even the lowest perturbative 
order for the velocity, equation (j67l) is a very complicated 
expression. We know from the numerics that, in the sta- 
tionary state, the velocity is given by a periodic function 



1 


1 ' 1 ' 1 ' 1 ' 


- 


perturbative solution (a) ■ 




full numerics (a) 




perturbative solution (b) 


1 


full numerics (b) 
1.1,1,1. 







FIG. 5: Comparison of the full numerical solution of equation 
(|47|l with the result obtained from the first non-vanishing per- 
turbative order eq. (|67p for (a) h = 3.0, c = 3.0, r] — 1.5 and 
(b) h = 1.0, c = 1.0, rj — 0.6. The units of t and z are chosen 
such that to = £ = 1. 



with periodicity u!~^. This recommends to aim a har- 
monic expansion of the mean velocity v, i.e. to ask for 
the Fourier coefficients ajv and in the ansatz 

00 

v{t) = ^ [gn cos Nut + bN sin Nu}t], (68) 

N=l 

It is now an important question, whether the first-order 
perturbative correction can further be simplified to get an 
analytic description of interesting features of the trajec- 
tory. One idea could be, to take only the lowest Fourier 
modes. In this section, we want to investigate, under 
which circumstances this could be possible. 

Starting from the first order result for v, given by eq. 
((67| . we express the disorder correlator by its Fourier 
transform 



A'[Zit)-Z{t')]= I g(ig)A(g)e'^^ 



[sin ujt— sin ujt'] 



(69) 



and expand the exponential term in a double Fourier se- 
ries in t and t', respectively: 



Aa sin ut 



n— — 00 
00 

^^,^_c(t-t')-iasi„c.t' ^ J„(-a)- 



c -I- inui 



Here, J,i(a) are the Bessel functions of the first kind. As 
we are interested only in the behaviour for large enough 
times (the steady state solution), we remove all terms 
that are damped out exponentially for t ^ from the 
very beginning. Note, that is indeed the time scale 
for the transience, as has been claimed before. 
For the mean velocity, we obtain 
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v{t) — hcosivt + if 



E /g(i<z)A(9)J„ 



q— Jn 



h\ e'(™+")"*(c-mw) 



(70) 



In principle, this is already a Fourier series representa- 
tion, not very elegant, though. The argument (to + n)ujt 
of the expansion basis exponentials promises a rather 
complicated structure for the coefficients. A first observa- 
tion, however, can already be made: Under the q integral 
we find an odd function {iq)A{q) and a product of two 
Bessel functions of order m and n, respectively. For the 
q-integral to result in a finite value, a function is required 
that is not odd in q. This necessitates the product of the 
two Bessel functions to be odd, or, equivalently, m + n 
to be an odd number. Whence, we conclude, that to first 
perturbative order, our symmetry argument (Fourier co- 
efficients for even N must vanish) is fulfilled exactly. 

It requires some tedious algebra to collect all contri- 
butions belonging to a certain harmonic order from the 
double series. Eventually, we obtain a series expansion 



/i = 2.7; c= 1.0; ?; = 0.6 



vit) 



cosojt - 



E 

N=l 



A 



N 



rj uj 



uj£ 



cos NLdt+ 



Sat I -, -, — I sinNujt 
h c Lot 



(71) 



Note, that taking w ^ is forbidden here, as we used 
cj 7^ while deriving the coefficients and moreover per- 
turbation theory breaks down (recall that h/rj ^ c/w). 
The same holds for £ — > 0. The remaining extreme limits 
cj — > cx) and i ^ oo are not interesting, since in these 
limits the disorder is rendered unimportant. Therefore, 
in the following, we assume finite (positive) values for £ 
and Ld and moreover set them equal to one uj — I ^ l,hy 
appropriately choosing the units for z and t. 

Now, we are left with three dimensionless parameters: 
h, c and rj. The dependence of the first order perturbative 
Fourier coefficients on 77 is trivial. The dependence on c 
is also evident, as can be read off from ([70)1 . For larger c, 
the system is more tightly bound to the non-disordered 
solution, supressing perturbative corrections. 

The most interesting but also the most difficult is the 
dependence of the Fourier coefficients on h. Actually, 
there are two competing effects. On the one hand, large 
driving strengths render the disorder unimportant in all 
cases accessible through perturbative methods. On the 
other hand, if one thinks of g{Z{t)) as a function of time, 
the more rapid Z{t) changes the more g fluctuates on 
short time scales and thus brings higher frequency con- 
tributions to v{t). The first remark is reflected in the 
overall weight of the Fourier coefficients as corrections 
to the non-disordered case, decreasing with h. The sec- 
ond idea is expected to express itself in the decay of the 
Fourier coefficients with N . The larger h, the weaker we 
expect this decay to be. 




FIG. 6: Plotting the logarithms of |^jv| and \Bm\ reveals the 
exponential decay with A''. In the regime where numerical 
errors do not dominate the result, a linear regression seems 
appropriate. 



c = 1.0 ; ri = 0.6 




FIG. 7: Performing the linear regression for many h yields 
slopes a and /3 appearing to depend on /i in a power-law 
fashion. 



The dependence of the higher harmonics on h is hid- 
den in the Bessel functions in (|70p . The first extremum 
of the Bessel functions shifts to larger values as to and 
n increase, respectively. However, the complicated way 
in which these Bessel functions enter An and Bn hin- 
ders an analytic access to the decay law. A numerical 
determination of the Fourier coefficients for the pertur- 
bative result reveals an exponential decay, as shown in 
fig. [HI The noisy behaviour for iV > 40 is due to nu- 
merical fluctuations. Note, that these fluctuations are of 
the order 10~^^, which is quite reasonable. The plot in 
fig. [S]is mere illustration of a more general phenomenon. 
This exponential decay has been found for many sets of 
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TABLE I: Results for the regression (f73)) . 



c 


V 


Cq. 


Cl3 




^0 


1.0 


0.6 


1.52 


1.52 


0.61 


0.61 


1.5 


0.6 


1.56 


1.56 


0.58 


0.59 


2.0 


0.6 


1.56 


1.59 


0.58 


0.60 


2.5 


0.6 


1.63 


1.62 


0.61 


0.61 


3.0 


1.0 


1.66 


1.66 


0.63 


0.63 


o.O 


l.U 


1 71 
1.(1 


l.Do 


U.DZ 


U.DZ 


4.0 


1.0 


1.69 


1.65 


0.61 


0.60 


4.5 


1.0 


1.72 


1.68 


0.62 


0.62 


5.0 


2.0 


1.71 


1.67 


0.61 


0.60 


5.5 


2.0 


1.73 


1.73 


0.61 


0.62 


6.0 


2.0 


1.77 


1.72 


0.62 


0.62 


6.5 


3.0 


1.76 


1.77 


0.62 


0.62 




c 



FIG. 8: Plot of the change of the prefactor C(c) = (Cc+C/3)/2 
in (|73|) on c. The linear fit yields a fairly tiny slope. 

parameters, thus one is led to the ansatz 

|A^|~|Je-"^ ; |i?A.|^^e-^^ (72) 

where a and l3 can be estimated through a linear regres- 
sion up to a suitable iV^ax- Of course, it is not expected, 
that a and /3 are distinct, nor that they depend on the 
parameters in different ways. Determining both just dou- 
bles the amount of available data. 

As our results are first-order perturbative, a and (3 
must not depend on rj. The main interest now focusses 
on the dependence of the decay constants on h. The 
results from a linear regression for a series of /i- values, c 
and 7] kept fixed, suggest a power-law dependence 

a{h,c) ^Co.{c) ■ h-^", l3(h,c)^Cp(c)-h-^^. (73) 

Fig. [7] displays this relation for a particular example. 
Repeating this data collection and subsequent regression 
for different values for c and 77 yields the results sum- 
marised in table ID While the exponent ^ appears con- 
stant ^ ft! 0.6, the prefactor seems to depend on c. An 
attempt to redo the same procedure, done for h, with 
the parameter c to gain information about the functional 



dependence of a and /? on c yields a complicated but 
rather weak dependence, which gives no further insight. 
The linear fit in fig. |8] gives a fairly tiny slope, so the 
dependence of the decay constants on c may be assumed 
to be weak. 

Certainly, it is desirable to ascertain the validity of this 
decay law beyond perturbation theory. In a few words, 
it ought to be explained, why we have not been able to 
do it. First of all, the logarithmic plots of the Fourier co- 
efficients in fig. [6] exhibit fluctuations around the linear 
decrease. This "noise" is authentic and not attributed 
to numerical inaccuracies. The exponential decay of the 
Fourier coefficients is superimposed on a true, compli- 
cated dependence. Hence, it requires a lot of data points 
to obtain reasonable data. Since the Fourier coefficients 
for even N vanish, in the example of fig. [5] the regression 
can be carried out over around 15-20 data points. This is 
a fair number. The quality relies heavily on the accuracy 
of the numerical determination of the Fourier coefficients. 
In Fourier analyses of the numerics for the full equation 
of motion (|47)) . we did not manage to get a precision bet- 
ter than of the order of 10~^. This means, the regression 
has to be stopped at iVmax, where log^Ar^^^^^ sa —7. In 
the example of fig. [SI this leaves us with less than 5 
data points. In view of the natural fluctuations, a linear 
regression is not sensible any more. 

In summary, we have numerically established the de- 
pendence of the decay constans for the Fourier coeffi- 
cients on h: a, 13(h) = C(c) • with ^ ~ 0.6. 



IV. CONCLUSIONS 

We have seen that the perturbation expansion for ac- 
driven elastic interfaces in random media fails for inter- 
face dimensions D < A. We have resolved this puzzle, 
and the reason for the strange behaviour of perturbation 
theory has been found to be connected with the disorder- 
dependent time scale r (cf. ((29)) ) which measures the 
time for the initial roughening of interfaces in random 
environments. Due to its appearence within the disor- 
der averaged velocity w as a transience relaxation time, 
attemps to determine v via a perturbation expansion in 
the strength of the disorder entail terms of unbounded 
growth in time. Although, most probably, the expan- 
sion yields the true solution if it could be summed up, it 
is useless, because the finite perturbative orders grow in 
powers of time. They offer a good approximation to the 
full solution only on time scales small compared to the 
transience relaxation time and fail to describe correctly 
the behaviour on large time scales. On the other hand, as 
we have signified, the perturbation expansion for D > 4 
is regular. 

Therefore, theoretical work towards an understanding 
of the important problem of ac-driven interfaces exposed 
to disorder has to take the route via the related mean- 
field problem. The mean-field equation of motion for- 
mally corresponds to = c» and admits a regular per- 
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turbative treatment that, where apphcable, agrees very 
weU with the numerics for the full equation of motion. It 
has been shown, that non-regular diagrammatic contri- 
butions cancel among each other, leaving a well-behaved 
perturbative expansion. 

Further, the solutions to the mean-field equation of 
motion have many features in common with the prob- 
lem in finite dimension. Therefore, they can be useful to 
study some properties of the original problem, like the 
velocity- hysteretis. 

The mean-field perturbation expansion helped to im- 
prove numerical results, which allowed us to establish the 
dependence of the decay constants of the Fourier modes 
on ft, as a power law. 



The three contributions to the action for this functional 
integral formula read 



So[z,z] = if z^,t[dt-rvl]z,,t (A7) 

Jx,t 

Sh[z,z] = i [ z^^tKt) (A8) 

Jx,t 

Sdis[z,z]^ — / Zx^tZx,t'MZx,t - Zx,t']- (A9) 

^ Jx,t,t' 
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Appendix A: The functional integral approach 

Apart from a perturbative expansion of the equation 
of motion ([2]), it is also possible to compute the disorder 
averaged correlation functions via a functional integral 
approach jUli^. Starting from © 

dtZx^t = rV|zj;,t + h{t) + u ■ g{x, z^^t) = F[zrc^t] (Al) 

(with Zx,t being a short-hand notation for z{x^t)), we 
note that 

for any fimctional Writing 
5{dtZx.t - F[zx^t]^ 



Vz^^t exp 



ZxA(^tZx,t - F[Zx^t]) 



(A3) 



and using the usual cumulant expansion for Gaufiian dis- 
order, we find 



{A[zx,t]) = 4 / '^^x,t'Dzx^t Azx^t]e' 



S[z,z] 



(A4) 



S[z,z] = So[z,z] +Sh[z,z] +Sdis[z,z]. (A6) 



Certainly, we can also compute correlation functions 
A[zx,t, Zx',t']- For example, we re-obtain the propagator 
(HJ by calculating the response function 



■ {Zx,tZx',t')s^ = G{x - x' ,t- t'). 



(AlO) 



To set up a perturbative expansion in u, we decompose 
in the same way as in section III Bl z^^t = Zt + C,x,t with 
Zt — (h/uj) sin ujt. Instead of Zx,t we now have to consider 
the functional integral over (^^.t- Now 



S[zX] = Sq[zX] +Sdis[z,C] 



(All) 



/ Zx.tZx,t'^[Zt + Cx,t — Zf — Cx,t']- 

Jx.t,t' 

(A12) 



To compute correlation functions perturbatively in pow- 
ers of u requires to expand the normalisation Z as well. 
A small reflection shows that the lowest order contribu- 
tion coming from Z is of order 0{u^). Thus, if we want 
to calculate the velocity to order 0{u'^) we can ignore the 
M-dependence of Z and write 



V{t) = h{t) + dt {Cx,t) 

= h{t) + dt iCx.tSdisiz, C])s, + (^("')- (A13) 



Of course, we may only retain iSdis[5,C] up to terms of 
order 0{u^), i.e. 
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5dis[5,C]=-v / z,,tZ:,,f[^[Zt- Zt.\ + /S.'[Zt- Zt.\{C,,^t-(:x,t')\+0{u^). (A14) 

^ Jx,t,t' 

Averages with respect to Sq are GauBian, thus Wick's theorem applies and with (jAlOp we recover the result from 
equation (fT3|) 

v{t) = h{t) - ^dt I ^'[Zt^ - Zt^] [{C,xaZx'm)s,, {CxmZx'm)so ~ {Cx,tZx',t2)sa {CxmZx'm)so\ 

^ Jx',ti,t2 

= hit) - u^dt f A'[Zt, - Zt,]{-iGix - x', t - h)}{-iGix ~ x', h - t2)]) + 0[u^) 

J x' ,ti ,t2 

= h{t)+u^ l\'[Zt - Zr]G{Q,t - t') + 0{u'^). (A15) 



It is clear, that we can continue equation (jA13|) to higher 
orders in u by taking into account higher orders in u of 
exp(5dis) as well as higher order corrections from Z. This 
way, it is possible to rearrive at the graphical expansion 
presented in section UlCi albeit along a little more com- 
plicated route. 

Appendix B: The regularity of all perturbative 
orders in case D > 4 

Let us start with an example graph at which we demon- 
strate the steps, that are then generalised further down. 



We consider the following fourth order contribution to 
the correction of the disorder averaged velocity 



' ' = ^ 

G{x — xi,t — ti)G{xi — X2,ti — t2)G{x2 — X3, t2 — h)x 

5{x - X2)5{xi - a;3)(-A"[Z - Z2])A'[Zi - Z3] (Bl) 

t tl t2 



= j dh Jdt2 J dt3{-A"[Z - Z2])A'[Zi- Z3]x 



000 



d^x' G{x -x',t- tx)G{x' -x,ti^ t2)G{x - x' , t2 - h) (B2) 
x'G{x ^ x' ,t — ti)G{x' — x,ti ~ t2)G{x — x\ t2 — t^) = 

d^fcld^fc2 ^_r[fc^(t-tl) + (fci+fc2)"(tl-t2) + fc"(t2-t3)1 ^p.O^ 

r 



We can now assign momenta to the branches of the tree and we can read off from (|B3p . that 

^ka ■ «, (B4) ka = ki , ki3 = ki + k2 , k^ = fe. (B5) 
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Obviously, ki describes the momentum "flowing" from 
the root to vertex 2 and fc2 is the momentum "flow" from 
vertex 1 to vertex 3. Estimating the disorder correlator 
derivatives in equation (jB2[) by constants (that we can 
drop since they do not influence the asymptotics t — >■ oo), 
all that has to be done is the integration of the right-hand 
side of (jBSp over the 3 time variables ti, t2 and t^. We 
start with the outermost leaf, ^3, and have 



This is certainly finite for D > A. 

Now, we turn to the general argument for all trees 
of general (even) order 2q. Let Tq denote the set of all 
rooted trees with q vertices, and 'P{T) all possible un- 
ordered pairings of vertices of T S Tq. Let moreover Bt 
be the set of all branches of the tree T. We want to agree 
that a branch b — (61, 62) has 61 always closer to the root. 
Then, we have in general for every order 



E2 

/ 



1 _ e-rklt2 



Proceeding in the same manner for the remaining time 
variables, we finally find 



< 



d^fcld^fc2 



1 



(27r)2C r3fc2(fci +fc2)2fc2' 

(B7) 



{V2q) 



E 



Ter2, Pev{T) 



A 



T,P- 



(B8) 



Here, At^p is a single diagram, namely a tree T where 
all pairs (^1,^2) S P are connected by dashed lines for 
the Gaufiian disorder average. The disorder correlators 
that enter At^p can be estimated by constants | A^™^ 
Z"]\ < Cm^'^"^, that we drop in the following. They are 
finite and do not influence the behaviour of At^p as t ~> 
00. Denoting the root vertex by r, we estimate 



tt, 



(B9) 



u JM^J di^o ( n / d'^-p.) exp -ty: kuu, - 

b^BT n p^-f h^Bx 



+ i Xpi ('^bi .Pi + '^bl ,P2 ~ ^f)2 ,Pl ~ '^ti2 ,P2 ) 

peP beBT 

The integration over the remaining x-coordinates brings delta-functions for the momenta: 

t, 

/ „ r H^K 

\Atp\ < 



(BIO) 



_ 1 

/d k f ^ 
j^—^ j dtb, j exp - r ^ kl{tb, - tb, 

n (2'r)^<5f ((^fci.pi + '^6i,p2 - <5b2,pi - Sb2,P2) j 



(Bll) 



■ beBT 



The g — 1 delta- functions for the momenta mean, that the 
net out-flow (away from the root) of momentum from a 
vertex pi equals the net in-flow of momentum for the 
Gaufiian partner vertex p2 (i.e. (^1,^2) G P)- There 
is no delta-function ensuring this for the root and its 
partner vertex. However, this is not needed. The root 
itself has only out-flow of momentum and that has to 
be absorbed by its partner vertex to ensure the balance 
for all other pairs. This explains, why no exponential 
function involving x appears any more in (IbTTI) . As a 
result, we can assign a momentum hp to each pair P 3 
p = (pi,P2)) describing the net momentum transfer from 
the vertex pi to p2- This insight advises to choose the 



momentum associated to the q Gaufiian pairs as the q 
integration variables kp. The momentum assigned to a 
bond k}, denotes the total flow of momentum through this 
bond, determined by the source or drain properties of the 
bond boundary vertices bi and 62. The rest is now easy, 
if we follow the estimate for the time integrals from the 
equation (jB6l) in the example before. We end up with 



< 



n 

peP 



d^fcp 
(2^ 



n Tki 

beBr ° 



(B12) 



One first puzzle related to (jB12l) is obvious: there may be 
bonds, that carry momentum k^ = 0. This exactly hap- 
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pens for lines that connect to a subtree which has only With regard to their topological structure, these are the 
internal Gaufiian pairings, i.e. the whole momentum flow only diagrams that have problems in D > 4, other such 
remains inside the subtree. Such lines are problematic in diagrams are constructed from them by a different choice 
all dimensions, even in the mean-field theory, where they of the root. Incidentally, I?| and D| are topologically 
correspond to the curly lines in the graphical expansion equivalent as well. We can easily see that in each dia- 
for the mean- field equation of motion (cf. section UlI C|) . gram Df, . . . , Df, there are 3 bonds that carry the mo- 
The resolution of this problem by cancellations among mentum ki belonging to the pair of the root and the out- 
several such diagrams is technically a little easier for the ermost vertex. The momenta associated to the other two 
mean-field case, where we performed it (cf. section llll C 21 Gaufiian pairs, ^2 and fca, occur only within exactly one 
and appendix lEl) . The idea of the mean-field proof ex- bond. Thus, according to equation (|B12p the integration 
actly applies here as well, just in mean- field we do not over ki produces problems in D < 6, whereas integration 
carry the load of momentum integrals. The adaption of over k2 and ks is harmless in Z) > 2. One can now check, 
the proof itself is straightforward. that the sum ADf + 2D2 as well as the sum of 2D^ + D\ 
Another problem in connection with (jB12|) is the fact, is regular in Z? > 4. The factors that I have put count 
that we have only q momenta to integrate over but 2q—\ the incidence of each diagram in the expansion, 
bonds. So for every given order it happens for a couple of The mechanism of combining diagrams to improve the 
diagrams, that some momentum, kp say, appears to the properties with respect to integration over some ki and 
power k~^'^ in (jB12p . Thus, the integration over kp has make it worse for some other kj can already be demon- 
infrared problems in < 2q, but integration over the strated at graphs of order 4. If, from and Z?! the root 
other momenta works in Z) > 2. Actually, this problem and the upper vertex closest to the root are removed and 
can again be resolved by summing several such diagrams the outermost leaf connected to the new root, we have 
in multiple steps. After the first step, one obtains a re- the two diagrams of fourth order (cf. (|23p) 
suit, that has kp divergent in Z) < 2q — 2 but another 
momentum, kpi say, divergent in Z? < 4 instead of Z? < 2. 
So, viewed as a function of kp the intermediate result is 
better behaved, but worse for kpi . Let us take a closer 
look at this. We start with a glimpse on this issue for 
the expansion of (ug) i where we have checked that every- 
thing remains bounded in D > A. The expansion of (ug) 
contains, among others, the following diagrams 




\ 



(B13) 




(B14) 




(B15) 
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o 1.1 

f /"dtie-^'='(*-*i) /dt2e-"'(*i-*^)A"[Zi-Z2] /'dt3e-i^'=i(*i-*^)A'[Z-Z3] (B17) 

J ki ,k2 ^ ^ J 



/ / dti e-"?(*-*i) / di2 e-r'='(*i-*^)A"[Zi - Z2] / e-r'=?(*^-*^)A'[Z - Z3 

J ki ,k2 •J ^ 





(B18) 
(B19) 



5 = / / dii e-"?(*-*i) /" di2 e-i^'='(*i-*^)A'[Z - Z2] f At e-r'^'^*^-^) A"[Zi - Z{t) 

JkiM J J J 



t ii t2 

- / / dti e-^'^?^*-*^' / dt2 e-r'='(*i-*=)rA;? / dtae'^'^*^-*^) A'[Z - Zajx 



t2 

y dre-r'=i(*^-")A"[Zi-Z(T)] 



(B20) 



The modification of D\ relies on integration by parts for 
the integral over ^2- Upon inspection of the integrals, one 
sees, that the integral over k\ in the diagrams D\ and D\ 
is well-behaved for i — ?> 00 in Z3 > 4, while the integral 
over k2 requires only D > 2. The sum S = 2Df + Z?! 
on the other hand behaves vice versa. Of course, in ex- 
pressions of order 4 not much is gained with that, but in 
higher orders this procedure can be used to understand, 
why the sum of all diagrams of a certain order 2q involves 
g — 1 integrals over momenta that are finite as t — )■ 00 
in Z? > 4 and one that works in D > 2. Note, that in 
equations (jB17IIB20| it is important to retain the disorder 
correlators, because cancellations have to be exact. 

It is now easy to extend the calculations provided in 
(|BT7llB20)) to see that 4:Df+2D^ (cf. equations ([BTI and 
(|B14[) ) is a regular expression in I? > 4. The sum 4I?f + 
2D2 has the critical dimension for ki, the momentum 
associated to the Gaufiian pair made of the root and the 
outermost leaf, reduced so that it is regular in > 4 
at the price that now k^ also needs D > A instead of 
D > 2. For reasons of space, as the explicit expressions 
for the diagrams are rather huge, we do not provide the 
full calculation here. 

The more general idea about the cancellation among 
trees goes as follows. Because too many details have 
to be ascertained which tends to result in a huge load 
of technicalities, we will merely sketch the procedure by 
mentioning all intermediate steps without proving the 
implicit claims. First of all, we note, that any tree T 



involves at least one distinguished momentum, we denote 
it by kT, the integration over which is regular in £) > 2. 
In (|B17IIB20)) we have seen that it is possible for any tree 
Tq of any given order q, that is regular in > 4, to find 
partner trees T^,. . . ,T™ of the same order, all regular 
in _D > 4, such that their sum is an expression in which 
the momentum associated to the root has the property 
of kr- Thus, we will without loss of generality assume 
that for any tree that is regular in Z? > 4, kr is the 
momentum fiowing out of the root. All our arguments in 
the following will also hold for trees for which this is not 
the case, if we repeat it for all partner trees and take the 
sum. 

Let, for some tree T the integration over some mo- 
mentum kp associated to a Gaufiian pair p be problem- 
atic in D > 2No + 2. This can happen if along the line 
from pi to p2 there are iVo vertices which form the root 
of independent subtrees {Si} or if the momentum flow 
along the route from pi to p2 is interrupted by Nq inde- 
pendent (w.r.t. the disorder average) inner subtrees {Si} 
and continues at the Gaufiian partner of the root of those 
subtrees. The first case corresponds to the diagram Df 
given by (jB13l) and the second scenario is exemplified by 
Dl which is depicted in ((BT6t (both for iVo = 2). Of 
course, also a mixture of both events, Nq in total, has 
the same effect, as is illustrated by Z?2 and Z)|. If, as 
a kind of induction hypothesis, we assume that all inde- 
pendent subtrees Si, . . . ^ Sn^ are regular in Z? > 4, we 
can describe the scheme how to find all trees that have 



19 



to be added to T to yield an expression that is regular 
in _D > 4. Let us consider first Nq — 2. Let si be the 
root of 5*1 and s'l be its Gaufiian partner vertex in 5*1. 
As mentioned above, the integration over the momentum 
/cs, associated to [si.s'i), is regular in Z3 > 2. Without 
loss of generality, we assume that in T, si is a vertex 
on the non- interrupted path from pi to p2- Then there 
is another tree T', for which the flow of kp along the 
connection between pi and p2 is interrupted by the con- 
nection between si and s'^. The sum of T and T' is then 
regular in D > 4. The integral over kp in the sum T + T' 
has decreased its critical dimension by 2 at the cost, that 
now the integration over k^ needs 13 > 4 to be bounded 
for large t. 

So far, this is the idea, how the cancellation among 
trees works in 13 > 4 to give a regular expression. For a 
thorough proof, we would have to give evidence for each 
single intermediate step. After all, the practical bene- 
fit of such a detailed proof is little and no further insight 
can be expected. Thus, although a rigorous proof has not 
been established, a consistent picture of the behaviour of 
the perturbation expansion has emerged. In D > 4 all 
perturbative orders are regular, and in D < 4 our some- 
what crude estimate of the disorder correlator derivatives 
by constants gives bad results. 



where we have introduced 



t 

LrnA<l,t) =Jra Jn / dtidt^ e'-(™*^+"*=) 



d^fc 



{2n) 



D 



-Tk^{2t-t2-t2) 



d^fc 



Tk'^ + \nuj 



(C4) 



The behaviour for t — ^ oo is dominated by the leading 
term LQfi{q,t), which is given by 



2 / 



->D 



a; / r2(27r)^ 
2 / ih\ f 



dkk 



D-5 



1-e" 



1 2 



U! / 1 2 



(C5) 



The function A]j{x) has already been introduced in equa- 
tion (I5T|) . The sub- leading terms are given by 



Appendix C: Analysis of the bush graphs 



To analyse the term 

Ti = 1 dtidt2 A[Zi - Z2] I 



-rk^{2t-ti-t2) 



(CI) 



from equation (j26l) . we start with the decomposition of 
the function A[Z{ti) — Z{t2)] in a double Fourier series in 
ti and t2, respectively. Recall, that Z{t) — {h/u))smujt. 



A[Z{ti)-Z{t2)]^ [ A{Q)J2'^rn{qh/L0)Jn{-qh/^ 



LU)X 



(C2) 



Here, A is the Fourier transform of A and J™ is the 
Bessel function of the first kind. For symmetry reasons, 
only terms with an even value of m -I- n contribute. This 
gives 

Ti^ I A{q) Lofl{q,t) + L2mfl{q.t) + 



(C3) 



L2m,0 ^Jq ( — ) J2m ( — I X 



Thus, 

YL2„ifi{q,t) = 



1-e 



-rk't 



rfc2 



rk^ + i2muj 



(C6) 



Here, the function a'^{x,y) is given by 



{t/^,u;t). (C7) 



1-e" 



2 cos 2m-;/ -I- 2m sin 2m-;; — ^e p 
V f ^ V 

(p2/y2) _,. 4^2 



(C8) 



which, for x 00 behaves in the same way, as ao (x) (cf. 
dTSH . independent of m 7^ 0. 

The last term in (|C3I) does not require further consid- 
eration. For m, n 7^ 0, the function L„i n remains finite as 
t — >■ 00. There is no infrared problem with the fc-integral 
in equation (jC4[) any more. Recalling the two time scales 
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T and that we have encountered before (cf. (|29)) and 
([50)1 ) we thus have 



4-D Q 

t\ ^ uA- , ft 



u^A^ / 1 



Hereby, we have introduced 
t 



2 I 9^ 



(C9) 



(CIO) 



kn \ -,ujt ] = { J 



A(,) Jo ( ^) J.™ f ^ 



The second factor in (|26|) 
T. = |d^'A(^-^)[Z-Z']|^ 



(Cll) 
(C12) 



-rfc^(t-t') 



(C13) 



can be treated in the same way, hke the second order 
graph in section Hi B| by spUtting off the Fourier-O-mode 

A(^-i)[z(i) - zit')] = (C14) 

FoUowing the calculations in section III Bl with Fq {ojt) 
replaced by F^^^ ^\u!t), we arrive at 



,2£2(p-l) 



mA ? 










Pd 





Id ( 



(C15) 



where 



Id I 



S 



D 



aD{t/d)F^^''-^\ujt) (C16) 



(2^) 
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dp p 



(C17) 



The integral over p in is certainly convergent for any 
D in the limit t oo, and the integral over t' converges 
for any D > 0, since p{t,t — t') is a bounded oscillation 
around zero (without zero Fourier mode) in t' . 



Appendix D: The width of ac-driven interfaces 

Apart from the velocity of the mean position of an in- 
terface in a random potential, there is another interesting 
quantity that deserves investigation: the mean square 
deviation of a given realisation from the mean. More 
precisely, we refer to the quantity 

w^{{{z)-zf). (Dl) 

In the first order of the perturbation expansion, w reads 



w={{{Z + uCi) 



(D2) 



In the case of infinitely extended interfaces, this quantity 
measures thus the typical width of the interface. 

So, for infinitely extended domain walls, the typical 
width to first order in perturbation theory is given by 

(cf. m) 



{Cl) {x,t) = J dhdt2 A[Zih) - Z(i2)]x 



d^fc 



(2^)^ 



rfc-^(2t-ti-t2) 



(D3) 



Comparing this to Ti, given by equation (j27l) . we find 
(^^) = Ti. Using equation (|28l) . we thus have for the 
asymptotics t — )■ oo 



(D4) 



The function Ad, given by pip , remains bounded for 
D < A and grows logarithmically in its argument in case 
D = A. Thus, the growth of the perturbative estimate 
of w in time is given by the prefactor and logi for 
D < A and D — A, respectively. 

So, in contrast to the first order perturbative result for 
the interface's velocity, which remains finite as t — >■ oo for 
D > 2, the width of the interface indicates the correct 
critical dimension D ^ A already to first order. 



Appendix E: Regularity of the mean-field 
perturbation expansion 



In section lTlI C 2l we have analysed, how the unbounded 
contributions, contained in the two diagrams that involve 
a curly line, mutually cancel in the second non-vanishing 
perturbative order. In this appendix, we are going to 
explain how this cancellation process generalises to all 
orders in perturbation theory. As before, for simplicity, 
we work with the diagrams for the disorder-averaged ve- 
locity, that arise by just removing the curly lines from 
the root of the diagrams for {() (cf. equation ([S^). In 
a velocity diagram contributing to the n-th order (recall. 
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that only for even n the corrections are non-zero), any 
curly line connects two trees of order p and q (both even) 
with the restriction p + q — n. Both trees appear in the 
expansion of lower orders, namely p and q, respectively. 
In the following, we sketch an inductive proof for the 
claim that the unbounded terms originating from trees 
with curly internal lines cancel among each other. 

Let us assume, that for order n we have achieved to 
ensure regularity. For every unbounded tree T, there 
is thus a set , . . . , T° of, let us call them cancelling 
trees, such that T + T^ + ...+T'' is a regular, bounded 
expression in time. As a starting point for the induction, 
take n = 4, where the validity of the claim has been 
verified in section UlI C 21 It is now the task to validate 
the regularity for order n-f 2. First of all, we consider the 
process of attaching the root of a regular tree S (with no 
internal curly line) of order s by a curly line to a vertex v 
of another regular tree R of order r — n + 2 — sto obtain 
a new irregular tree A of order n + 2. The vertex v must 
be connected to another vertex w G Rhy a dashed line, 
to carry out the Gaufiian disorder average. Without loss 
of generality, we assume that v is connected to w by a 
path that first makes a step towards the root. The rules 
for the diagrammatic expansion ensure, that there is a 
maximal regular subtree T C R, which contains v and w. 

Using partial integration, it is possible to move the 
vertex to which S is connected (via the curly line) to a 
neighbouring vertex in T. Thus, it is possible to move 
the connection vertex along the unique way (in T) from 

I 



The order of the derivative (i.e. the number of outgoing 
lines) of w and v are denoted by fi and i/, respectively. 
The time, at which the whole diagram is to be evaluated, 
is t, the time corresponding to the vertex to which v is 
connected is given by Ti, ti is thus the time associated 
to V and so on. The time of w is r. Thus, we see, that if 
w is not the vertex to which v is directly connected (then 
Tl 7^ r in general), the first expression after partial inte- 
gration cannot be a valid diagram: v has lost one order 
of derivative {v — 1 lines go out instead of v), but the 
derivative of the correlator A has not changed. A valid 
diagram however reappears, when the connection of S 



V to w. We are going to show, that once w is reached, 
we have obtained the cancelling tree which is unique. 
Diagrammatcially, the process of moving the connection 
vertex from v io w reads: 




Here, the blank circle represents S, the light grey circle 
stands for the subtree Ri of i?, to which v connects and 
the darkgrey shaded circle denotes trees which run out 
of V (summarised in the following as i?2)- Certainly, in 
general there may be dashed lines between the dark- and 
the lightgrey circle, which we have omitted as they are 
not relevant for the forthcoming discussion. The dotted 
line just serves as a joker - it is not important to specify 
how many lines go out of v. The last term D collects 
the left-over terms from the partial integration. Note, 
that, if it takes several steps to go from v io w, the inter- 
mediate expressions (in the partial integration) are not 
in accordance with our diagrammatic rules, because the 
order of derivative of the disorder correlators does not 
appear correctly (it remains the same but the graph has 
changed). Keeping this small peculiarity in mind, it is 
nevertheless instructive to think in diagrammatic terms. 

To illustrate the procedure, we take a look at the first 
step: 



(E2) 



(E3) 



has reached w. Then, v has lost an outgoing line, but 
w received one more and we indeed have achieved a can- 
celling tree: the factor (—1)'' remains, the true diagram, 
however, has (—1)'^^^. The signs are different, thus the 
two trees cancel. The left-over term from the partial in- 
tegration is again regular. This can be seen because all 
time integrals carry an exponential damping term. It is 
clear, that this is generally true for every partial integra- 
tion step. 

To go one step further, we assume now S to be irregu- 
lar. Essentially, the same procedure works, but there are 
more cancelling trees: one has take all cancelling trees 



'^^-^ =Ri{t) Jdti 




Tl Tl 

Ri{t) I At2S{t2) I dtie-^(^i-*i)(-l)'^A(^+'')[Z(T)-Z(ti)]i?2(ii) 



Tl tl 

i?i(t) j Ah e-^^'-'^^fSih) J dt2 e-^(*i-*^)(-l)''A(^+^)[Z(r)-Z(t2)]i?2(t2) 
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{S^} for S into account (which exist by induction hy- 
pothesis), thus S is replaced by and thence the 
left-over terms are again regular. 

A possible irregularity of R can be accounted for in 
the same way. It is, however, important to explain why 
this is possible, i.e. what are v and w in the cancelling 
trees for R. In the case of irregular S the problem was 
easy, since all trees have a unique root. As we have seen 
already, the procedure of creating cancelling trees does 
not change the structure of regular subtrees. Thence, all 
cancelling trees for R contain T. This makes clear, which 



V and w have to be chosen in the cancelling trees: they arc 
well-defined in T and T is a well-defined subtree of the 
cancelling trees. Thus, repeating the whole procedure 
described above for all cancelling trees of R yields the 
complete set of cancelling trees for A in the most general 
setting. 
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Symbol Quantity Reference 

z{x, t) interface profile Sec lII Al 

D internal interface dimension Sec lII Al 

h amplitude of the driving force Eq.Q and Eq.© 

CO frequency of the ac-driving force Eq.([l]) 

r elastic stiffness of the interface Eq.(15]) 

c elasticity constant (mean- field) Sec lIII Al and Eg. (1471) 

£ disorder correlation length Sec lII Al and Eq.(0 

u disorder strength Eq.Q 

?7 disorder strength (mean- field) Sec lIII Al and Eq.((2]) 

A disorder correlator Sec lII Al and Eq.((S]) and Sec lIII Al 

g disorder configuration Eq.((2| and Ea. (|47p 

Lp Larkin length Eq.© 

Z{t) solution in absence of disorder Sec lII Bl and Sec lIII CI 

disorder correction to the solution Sec lII Bl and Sec lIII Cl 

G propagator Eq.® and Ea. (157|) 

A inverse smallest length scale (UV-cutoff) Sec lII Bl 

r roughening time Sec lII Dl and Eq. (|29|I 

1? additional transience time scale Sec lII Dl and Eq. (|30|I 



